Finite Temperature Phase Diagram of a Two— Component Fermi Gas with Density 

Imbalance 
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We investigated possible superfluid phases at finite temperature in a two-component Fermi gas 
with density imbalance. In the frame of a general four-fermion interaction theory, we solved in the 
BCS region the gap equations for the pairing gap and pairing momentum under the restriction of 
fixed number densities, and analyzed the stability of different phases by calculating the superfiuid 
density and number susceptibilities. The homogeneous superfiuid is stable only at high tempera- 
ture and low number asymmetry, the inhomogeneous LOFF survives at low temperature and high 
number asymmetry, and in between them there exists another possible inhomogeneous phase, that 
of phase separation. The critical temperatures and the orders of the phase transitions among the 
superfiuid phases and normal phase are calculated analyticaUy and numerically. The phase diagram 
we obtained in the temperature and number asymmetry plane is quite different from the one in 
temperature and chemical potential difference plane for a system with fixed chemical potentials. 

PACS numbers: 03.75.Ss, 05.30.Fk, 74.20.Fg, 34.90. +q 



I. INTRODUCTION 

Recently, the superfluidity in a two-component Fermi 
gas with population imbalance promoted great interest in 
both experimental and theoretical studies. The key point 
is the ground state of fermion pairing between different 
species with mismatched Fermi surfaces, and it is directly 
related to the studies of atomic Fermi gas with population 
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16, 17, ISl, superconductivity with Zeeman splittin, 
20, l2l|, |22|, isospin asymmetric nuclear matter [23 
2M and dense QCD matter [H [13, [H, [H, H^. In the 
investigation, two exotic states, the homogeneous and 
isotropic Sarma (lOj or breached pairing [sil [32| state and 
the inhomogeneous and anisotropic Larkin-Ovchinnikov- 
Fudde-Ferrell (LOFF) [13, HU, [12 [H state, are especially 
concerned. While these two interesting states have never 
been observed exp er iment ally, the recent progress in the 
study of atomic Fermi gasjsj, [H, [s^ may provide us a 
way to realize and observe them. 

There are two crucial problems related to the Sarma 
state due to the existence of gapless fermionic exci- 
tation. One problem is its thermodynamical instabil- 
ity compared with the fully gapped BCS state, when 
the effective chemical potentials for the two species are 
fixed. This is called Sarma instability [l^. It is now 
accepted that the Sarma instability can be avoided in 
systems under some physical constraints such as fixed 
particle numbers [H, [s^ instead of fixed chemical po- 
tentials, long range gauge interaction with charge neu- 
trality requirement '2 7], a proper finite range attrac- 
tive interaction[32], and in the strong coupling BEC 
regionj^ [sl]. The second problem for the Sarma state is 
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its negative superfluid density [1,[3^ and negative number 
susceptibility 3 . The former leads to negative Meissner 
mass sauared|40i [Zlj for charged systems and indicates 
that the LOFF state is probably energetically favored 
than the Sarma stateji^Ii^) and the latter shows that 
the superfluid-normal phase separation iJj is also more 
stable than the Sarma state. 

The above discussed negative superfluid density and 
number susceptibility are obtained in weak coupling BCS 
region and at zero temperature. Recently, it is argued 
that the Sarma phase will be free from negative super- 
fluid density and number susceptibility in strong coupling 
BEC region [1, [1, [sl], and both experimental and theo- 
retical studies support that the normal-superfluid phase 
separation is energetically favored around the unitary 
region d, [3^. 

In this paper we are interested in the BCS region 
where the LOFF phase is shown to be energetically 
favoredj^ [To| . An important and interesting phe- 
nomenon in this region is the enhanced pairing corre- 
lation at finite temperature [H, [4^. The pairing gap 
is not a monotonous function of temperature and its 
maximum is located not at zero but at finite temper- 
ature. Especially, for a large number asymmetry, the 
superfluidity appears only at finite temperature. This 
strange and interesting phenomenon was found in the 
studies of isospin asymmetric nuclear matt er (23l [2^ , two 
flavor color-superconducting quark matter [27|, and gen- 
eral breached pairing superfluidity [sj- However, it is 
recently found that the superfluid density of the homo- 
geneous Sarma phase is negative at low temperature and 
becomes positive only in a temperature window near the 
critical temperature Tc[43, [4^. As a result, inhomoge- 
neous phases such as LOFF phase can enter the low 
temperature region, and the phase structure should be 
re investigated. 

The purpose of this paper is to analyze the phase struc- 
ture of a two-component Fermi gas with density imbal- 
ance in the BCS regime. Our analysis is based on the 
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mean field approximation which is believed to be a good 
treatment in the BCS regime even at finite temperature. 
We take LOFF phase into account. For the sake of sim- 
plicity, we consider the simplest pattern of LOFF phase, 
namely, the single plane wave LOFF state or the so-called 
FF state. While the consideration of a general LOFF 
state is complicated, we will argue that the phase struc- 
ture with the FF state will not be changed qualitatively. 
Due to the thermal excitation, we will not distinguish 
between the gapless and gapped phases and simply call 
them homogeneous superfluid (HS). 

After a brief review of the theoretical frame for the 
attractive two-component Fermi gas with population im- 
balance in Section ini we calculate in Section Hill the su- 
perfluid density and number susceptibility as a function 
of temperature at fixed population imbalance and obtain 
the turning temperature where the superfluid density and 
number susceptibility both change sign. In Section HVl wc 
include the LOFF state and separate the homogeneous 
region from the possible LOFF region. By calculating 
in Section |V] the number susceptibility x for the LOFF 
state, we further distinguish the stable LOFF region at 
high population imbalance from unstable LOFF region 
at low population imbalance. We obtain the phase dia- 
gram in the temperature and population imbalance plane 
in Section |Vl We summarize in Section IVTl The natural 
unit of c = /i = /cb = 1 is adopted through the paper. 



II. FORMALISM 

The physical system we are interested in in this pa- 
per is an infinite system composed of two species of 
fermions with attractive interaction in three dimensional 
free space. Generally, the system can be modelled by the 
Lagrangian density 



where ip„ (x) are fermion fields for the two species denoted 
by t and | with space-time a; = (t,x), g is the coupling 
constant, m is the fermion mass, and /i| and are the 
chemical potentials. 

For attractive coupling g we can perform an exact 
Stratonovich-Hubbard transformation to introduce the 
pair field <i> ffV'iV'T its complex conjugate <i>* ~ 
gxlAxp]. With the Nambu-Gorkov field defined as 5* = 



(^l,"!/"!)^, the partition function can be expressed as 



in the imaginary time (r = it) formalism of finite tem- 
perature field theory, where /3 is the inverse temperature, 
f3 = 1/T, and the kernel K. is defined as 
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Integrating out the fermionic degrees of freedom, we ob- 
tain 



with the effective action 



Seff[<i>*,<S>]-- / dr Syi Trln/C[$*,$]. 
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(4) 



(5) 



For a dilute gas, we can replace the bare coupling 
constant g by the low energy limit of the two-body T- 
matrixQ, 
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with the s-wave scattering length and fermion energy 
= P^/(2m). 

For an arbitrary scattering length at finite tempera- 
ture r, we should take into account the contribution from 
the pair fiuctuations and pseudogapjij to the thermo- 
dynamics of the system. Since in this paper we focus on 
the BCS region where the scattering length is negative 
and small and the coupling is relatively weak, the effect 
of pair fiuctuations and pseudogap can be approximately 
neglected, and the mean field approximation is a good 
approach to investigate the phase structure. In the mean 
field approximation, we replace the pair field $ and its 
complex conjugate by their expectation values. To have 
a unified treatment for both the homogeneous and LOFF 
superfiuid, the order parameter for the superfluid can be 
defined as 

{<^{x)) = Ae^^i-" , = Ae^^^'i-", (7) 

where A is the amplitude of the order parameter and 
can be taken to be real, and 2q is the pair momentum 
in single plane wave LOFF state. Obviously, q = and 
q 7^ correspond, respectively, to the HS and LOFF 
states. 

Since a general LOFF state can be considered as a 
superposition of single plane wave LOFF states, it corre- 
sponds to a deeper minimum of the free energy of the sys- 
tem, compared with the single plane wave LOFF state. 
Therefore, a general LOFF state should be more stable 
than the simplest LOFF state, and the simplest LOFF 
phase, if it exists, will be replaced by a general LOFF 
phase. In this sense, while the details of our phase dia- 
grams obtained in the following will be changed, the qual- 
itative phase structure of the system will remain, when 
a general LOFF state is included. 

After a phase transformation for the fermion flelds, 
Xa — e*'i"'0CTj the thermodynamic potential in mean fleld 
approximation can be evaluated as a summation of quasi- 
particles[43| 
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where Ea and Eb are the quasi-particle energies 



Ea 
En 



Sji 
6ji 



p • q/m 



(9) 



and 



with = (p2 + q2) /(2m) Ep ^ + 

average chemical potential /i = + /^x) ^'iid chemical 
potential difference 6fi = — /Xf ) /2. Note that, since 
the replacement of the bare coupling g by the s-wave 
scattering length Us , the ultraviolet divergence in the first 
momentum integration in ((8]) is removed, and there is no 
need to introduce a momentum cutoff. 

The thermodynamic potential obtained is a function 
of T, /i| and fii with A and q initially undetermined or- 
der parameter and pair momentum of the superfluid. In 
the spirit of thermodynamics, the physical system is de- 
scribed only by T, fi^ and /i| , and A and q as functions 
of T, /i| , are determined by the minimum thermody- 
namic potential. The calculation of the first order deriva- 
tives of O with respect to A and q gives the coupled gap 
equations. 
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ifiEs) - /{Ea)) 



E, 



(1 - /(Ea) - /{Eb)) 



,(10) 



where we have chosen a suitable frame with the z axis 
along the direction of the pair momentum, q = (0, 0, q), 
and f{x) is the Fermi-Dirac distribution function. It is 
easy to see that q = is a trivial solution of the gap 
equations, which corresponds to the homogeneous and 
isotropic phase. For a system with fixed numbers of 
species, the chemical potentials and in (fTO|) are 
solved from the fermion number densities n| and ni de- 
rived by the thermodynamic relations — —dVt/djia, 



Sp 

{2tt 



Mp), (11) 



where n-\{p) and ni{p) are the occupation numbers for 
the two species. 



"t(p) = ulf{EA)+vlf{-EB), 
niip) = ulf{EB)+vlf{-EA) 



(12) 



with the coherent coefficients Up = (l + ^p/Ep) /2 and 
vl = il-^p/Ep)/2. 

The gap equations (jlOp together with the fermion num- 
bers ([TT|) determine the physical values of A and q as 
functions of n-f and n|. Generally, there exist a homo- 
geneous solution with q = and a LOFF solution with 
q 7^ 0. To determine which one is favored at fixed fermion 



numbers, we should compare their free energies defined 
as 



(13) 



and the lower one corresponds to the ground state of the 
system. 

To explicitly describe the asymmetry between the two 
species, we use the total density n = nij + ni and the 
number asymmetry parameter a = (n^ — 't-|)/(?ix + n^) 
as variables instead of n-f and n^. Without loss of gener- 
ality, we suppose a > in the following. Like a free Fermi 
gas, we introduce a Fermi momentum pp or a, Fermi en- 
ergy Ep = p|,/(2to) through n = p|,/(37r^). If we scale 
all the variables with energy dimension by Ep and the 
variables with momentum dimension hy pp, the solution 
of the gap equations depends only on three dimensionless 
variables, the coupling pj^Cs, population imbalance a and 
the scaled temperature T/Ep. For the following numer- 
ical calculation in the BCS region with — 1 < ppUg < 
or — oo < \/{ppas) < —1, we take ppag — —0.6. This 
is a typical value for systems with BCS pairing. For 
example, for the pairing between ^Li atoms in states 
\F = 3/2, mF = 3/2) and \F = 3/2, mp = 1/2) with 
scattering length as = —2160aB where as is the Borh 
radius and typical density n = 3.8 x lO^^cm"'^, we have 
ppUs = —0.56. 



III. TEMPERATURE BEHAVIOR OF THE 
HOMOGENEOUS PHASE 

In this section we focus on the homogeneous phase. 
We analyze the temperature behavior of the pairing gap, 
the superfluid density and the number susceptibility, and 
show that there should be a phase transition from homo- 
geneous phase at higher temperature to some inhomoge- 
neous phase at low temperature. 



A. Order Parameter A 

At zero temperature, the gap and the chemical po- 
tential difference for an asymmetric system with a ^ 
satisfy the relation > A, which can be seen directly 
from the number equations. This means that the homo- 
geneous phase at T = must be a gapless superfluid [4l| 
with the gapless branch Eb- In the BCS region, we 
have a positive chemical potential fi which is not far 
from the Fermi energy Ep, and the chemical poten- 
tial mismatch 5fi is much less than fi, the two gapless 

nodes occur at momenta pi = \J 2m(/i 



A2) 



and p2 = \j2m[^i -f Sp,'^ — A^). The behavior of A is 
quite different from the one for symmetric BCS gap Aq 
with the same coupling ppUs- In weak coupling, there is 
an analytical relation between themjiol, 



A{Sfi) = y^\;{26j7^Ao), Ao/2<Sfi<Ao (14) 



4 



in terms of the chemical potential mismatch or [47 



B. Superfluid Density ps 



A(a) = Ao-\/l — a/aoQ{ao — a) 



(15) 



in terms of the population imbalance, where Q{x) is a 
step function and the critical imbalance ao is given by 
ao — 3Ao/(4i?i?). The solution Sji — Aq corresponds to 
a — Q and S^m = Ao/2 to a = ao. 

For the symmetric system, it is well-known that the 
pairing gap Aq will be suppressed by thermal motion at 
finite temperature, as shown by the solid line marked 
with a = in FiglT] However, this monotonous tem- 
perature behavior is not always true for the asymmetric 
system. The temperature effect not only deforms and re- 
duces the mismatched Fermi surfaces which leads to the 
usual suppression of the gap, but also makes the overlap 
region of the two species wider which favors the conden- 
sate. The competition of the two opposite effects results 
in a non-monotonous temperature behavior of the pair- 
ing gap. In Fig[Tl we show the gap A as a function of 
temperature for different values of population imbalance 
a and at the typical BCS coupling ppds = —0.6. In the 
case with weak imbalance, while the superfluidity starts 
still at T = 0, it firstly goes up with increasing temper- 
ature, then drops down, and finally reaches zero at the 
critical temperature Tc- In the case with strong imbal- 
ance, the superfluidity starts even at a finite temperature 
To > and vanishes at Tc- The phase transitions at the 
two critical temperatures To and Tc from normal phase 
to HS is of second order, and the critical behavior of the 
gap can be conventionally expressed as[48| 



A(r) = ipTc{l~T/Tcf\ 



T 
T 



T^ 

^ O 

t: 



(16) 



where To, Tc and the dimensionless coefficients (p and Lp all 
depend on the population imbalance a. When a exceeds 
the maximum value aus for the HS phase, there is no 
more homogeneous superfluidity at any temperature. In 
FigUwe have ans = 0.059. 



When 
locity 

($)e2' 



the superfluid moves with a uniform ve- 
Vs, the condensate transforms like ($) 
>v. x^ ($*) ^ ($*)g-2^mv,■x^ and the supercurrent 



is and the superfluid density Ps are defined via 



f2(v,) = r!(0) 



1 2 



(17) 



From the introduction ([7]) for the pair momentum q and 
its gap equation, the supercurrent keeps zero. 



is = 



dvs 



= 0, 



(18) 



and the superfluid density ps is directly related to the 
pair momentum susceptibility 



q=0 



(19) 



From the stability condition d^il/dq^ > for the pair 
momentum q, the superfluid density controls the stability 
of the homogenous phase. When ps is negative, the HS 
state is unstable. It is necessary to note that, up to the 
square term, the expansions of the free energy J-' in Vg 
and q are the same as the expansions of Jl. Therefore, we 
can safely use the deflnition (fT9l) to analyze the superfluid 
stability for the system with fixed number densities [49]. 

Using the result in [s^, [4l[ , the superfluid density at 
finite temperature can be explicitly expressed as 
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(20) 



with the definition f'{x) — df{x)/dx. Note that the inte- 
gration term is always negative, and the asymmetry be- 
tween the two species may induce a negative superfluid 
density ps- At zero temperature and in weak coupling 
limit, Ps can be approximately expressed as[4l|. 




FIG. 1: The pairing gap A, scaled by the Fermi energy Ef, 
as a function of T, scaled by the Fermi temperature Tp = 
Ep/kB, for different values of number asymmetry a. 
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5iiQ{5p. - A) 
^J5p? - A2 



(21) 



It is easy to see that in the Sarma phase with A < ps 
becomes negative. This is the so called magnetic insta- 
bility, since it is directly related to the negative Meiss- 
ner mass squared ra\ = e^ps/m^jJll if the fermions 
are charged. This dynamical instability implies that the 
LOFF state with finite pair momentum has lower free 
energy than the Sarma state at zero temperature. 

From the definition, the superfluid density vanishes at 
the critical temperature Tc for any imbalance < a < 
CiHS and also at the other critical temperature To > 
when the imbalance is close to the maximum ans- I^i 
weak coupling case, the chemical potentials and pi 
can be safely regarded as the Fermi energies of the two 
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species, and near the critical temperatures and To the 
superfiuid density behaviors as[48| 

= - T/TJ > 0, T^T-, 
Ps^Po{l-T/To)<Q, T^T+, (22) 

where pc and po depend only on the population imbal- 
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enough, ps is zero till the starting temperature To of the 
superfluidity, then becomes negative first and turns to be 
positive at T\. In Fig[31 we plot the two critical temper- 
atures To and T^ and the intermediate temperature T\ as 
functions of asymmetry parameter a. The temperature 
T\ separates the stable HS phase with > at high 
temperature from the unstable HS phase with pa < at 
low temperature. The negative superfiuid density indi- 
cates that the LOFF state may be energetically flavored 
at low temperature, and there may be a phase transition 
from LOFF phase to HS phase at T = Ti . We will come 
back to this point in the next section. 
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FIG. 3: The two critical temperatures To and Tc and the 
intermediate temperature T\ as functions of the asymmetry 
parameter a in the BCS regime with ppas = —0.6. ps > 
0, X > and ps < 0, x < indicate, respectively, the stable 
and unstable HS phases. 
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FIG. 2; The superfiuid density ps (dashed lines), scaled by 
Po = mn, and number susceptibility x (solid lines), scaled by 
Xo ~ n/Ep, as functions of T, scaled by Tf, for three values 
of number asymmetry in the BCS regime with ppcis ~ —0.6. 



C. Number Susceptibility x 

The general stability condition for a two-component 
system against changes in the densities of its com pon ents 
is described by the total free energy of the svstem|50l. l5l| , 
F — J (fix!F[n^{x)]. Considering its fluctuations induced 
by small changes 6na{x), the first-order variation SF 
vanishes automatically due to the number conservation, 
/ d'^x^rio- (x) = 0, and the second-order variation 5^F is 
given by the quadratic form 



S^F^Ud^^ E 



-drirrSnrr' ■ 



(23) 



ance a. Up to this point, we have shown that the super- 
fluid density is negative at low temperatures but should 
be positive near the critical temperature. In Fig[2] we 
show the superfluid density ps as a function of T for dif- 
ferent values of population imbalance a. For the symmet- 
ric system with a = 0, ps is always positive in the whole 
superfluidity region. For systems with weak imbalance, 
Ps is negative at low temperature, then changes sigh at 
an intermediate temperature Ti, and keeps to be posi- 
tive at high temperature. When the imbalance is strong 



Therefore, to achieve a stable homogeneous phase, the 
2x2 matrix d^T jdn^dnai should be positively deflnite, 
namely, it has only positive eigen values. This matrix is 
hard to calculate since it is not easy to express explic- 
itly the free energy as a function of densities. However, 
from the relation between the free energy Tin^) and 
thermodynamic potential f2(/io.), it is easy to check that 
the stability condition to have positively deflnite matrix 
J- / dn^rdn^i is equivalent to the condition to have pos- 
itively definite matrix ^d^fl/dfiadpa' ■ 
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For systems without mass difference between the 
two species, the condition to have positive eigen- 
values of —d'^V./dficrdna' can be reduced to the 
constraint 's'l that the imbalance number susceptibility 
X = —{d^fl/ddfi^)fi ~ {dSn/ dSji) ^ should be positive. 
For X < 0, the density difference 5n = — rt| increases 
with decreasing chemical potential difference (5/i, which 
is certainly unphysical and means instability of the su- 
perfluid against the phase separation (PS). Employing 
the gap equation which determines the condensate as a 
function of chemical potentials, A = we can ex- 

press the number susceptibility % as a direct and indirect 
parts, 

fd5n\ fd5n\ f dA\ 



From the expression 



Sn = 



(27r)3 

which leads to 
/ dSn\ 

fdSn\ 



(25) 



^[f'iEA) + f'{Es)], 
'^'P ^ [f{EB)~nEA)],{26) 



and the gap equation (dil/dA)^ Sp, ^ which results in 



d 
dSfi 



namely. 



dA 



9A2 



/ dA 



(dA\ _ / dSn 



\dA^ 



we have 



X 



\d5^jL J 



f-L.SfJ, 



0, 

(27) 
(28) 

(29) 



Since (dSn/dSfj,)^,/^ is always positive, the stability 
condition is controlled by the gap susceptibility k — 
[d^Q./ dA'^)^^s^ which determines if the solution of the 
gap equation is the minimum of the thermodynamic po- 
tential. 

We now consider the relation between the gap suscep- 
tibility and the superfluid density. It is easy to explicitly 
write down the gap susceptibility. 



l~ f{EA)^ fjEs) 

En 



JM^eI 

{r{EA) + f'iEB)) 



9 



(27r)3 



E„ 



(30) 



where the last line vanishes automatically due to the gap 
equation in the superfluid phase. In weak coupling limit, 
the number susceptibility and gap susceptibility in the 
superfluid phase at zero temperature can be evaluated as 



X 



3n 
2E^ 



mpp 



6^id{Sn - A) 



(31) 



which are negative at any imbalance a. On the other 
hand, applying partial integration to the number equa- 
tions, the superfluid density can be expressed as 



d^p 



A 



2 r 



(27r)3 3 E^p 



1 - /{Ea) - f{Et 



En 



inEA)+nEB)) 



(32) 



Comparing k and ps in the superfluid phase, the only 
difference is the factor of p^/S in the integrand function 
of Ps- In the BCS region, the integration is dominated by 
a narrow momentum window around the average Fermi 
momentum po ~ ^/2mp, and we have 



Ps ^ -jn. 



(33) 



From this proportional relation between k and ps and 
note that k is the denominator of the second term of 
(j29p . while the superfluid density and number suscep- 
tibility have different amplitudes, they change sign al- 
most at the same position. We have checked numerically 
that such an approximation holds well in the BCS region 
— oo < l/{ppas) < —1. When ps changes from Ps > 
to < at the turning temperature Ti, x changes from 
positive infinity to negative infinity approximately at the 
same temperature. Therefore, when ps is not too large, 
namely when the second term of ((29)) is guaranteed to 
dominate x, ps and x have almost the same positive and 
negative regions in the BCS limit. This temperature be- 
havior of superfluid density and number susceptibility 
is clearly shown in Fig[2] where we presented ps and x 
as functions of temperature. From the stability analy- 
sis at different population imbalance a, we plot in FigO 
the phase diagram of the homogeneous superfluid in the 
T — a plane for the coupling ppCLs = —0.6. The tempera- 
ture Ti (a) separates the stable HS phase at high temper- 
ature where x and p^ are positive from the unstable HS 
phase at low temperature where x and ps are negative. 
Because of the divergence of x at the turning point Ti, 
the possible phase transition from PS to HS must be of 
first order. 

From the above analytic and numerical calculations in 
the BCS region, the conditions to guarantee a stable ho- 
mogeneous superfluid described by the superfluid density, 
number susceptibility and gap susceptibility are equiv- 
alent. We can show that this is also true for unequal 
mass system. For instance, the breached pairing state 
proposed in [S^l which is the minimum of the thermody- 
namic potential should also be stable against the LOFF 
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state and PS state. Certainly, when the coupling is strong 
enough, the difference among these stable conditions will 
become remarkable. Due to the factor of p^/3, the su- 
perfluid density is easier to be positive in strong coupling 
region than the number susceptibility, this is the reason 
why the authors in [3] found that the positive superfluid 
density is a weaker condition than the positive number 
susceptibility. At finite temperature, we found that for a 
strong enough coupling, the turning temperature Ti for 
the superfluid density ps is remarkably different from the 
one for the number susceptibility x- Generally, the super- 
fluid density becomes positive at a lower temperature, 
namely, there exists a large region where ps is positive 
but X is negative. However, when the coupling is strong 
enough, the naive mean fleld treatment is not valid at 
flnite temperature and we should employ some proper 
flnite temperature theory for BCS-BEC crossover [44j. 

In the BCS region considered in this paper, from FigsO 
and [HI the LOFF and PS are both possible phases of the 
superfluid at low temperature T < Ti, the question left 
is to determine which of them is stable. 



IV. LOFF PHASE VS HOMOGENEOUS PHASE 

We now take the LOFF state into account. Firstly we 
consider the LOFF state with a < ans where the HS 
phase can survive. Before numerically solving the cou- 
pled gap equations (fTO|) for A and q together with the 
number equations for and 71 j^, we discuss ana- 
lytically the pair momentum in the temperature region 
close to the turning point Ti. The first order derivative 
of the thermodynamic potential with respect to the pair 
momentum can be written as 



dq 



qW{q) 



(34) 



which gives the gap equation for the pair momentum, 
qW{q) = 0. The homogeneous and isotropic state cor- 
responds to the trivial solution 9 = 0, and the solution 
corresponding to the LOFF state is given by W{q) = 0. 
From the second derivative. 



dq 



2 ^W{q)+qW'{q) 



(35) 



and the relation between the superfluid density and the 
pair momentum susceptibility demonstrated in Section 
we have 



which means 



W{0) = 



(36) 



(37) 



at Ti where p^ changes sign continuously. Providing the 
LOFF solution is unique (we have checked this numeri- 
cally at least in the BCS region), the LOFF momentum q 



must vanish at Ti, and in the temperature region below 
and close to Ti, q is very small and we have approxi- 
mately the relation between the free energies for LOFF 
state and HS state, 



^LOFF — J'HS 



2m? 



0{q') 



(38) 



Therefore, in the unstable HS region with ps < at T < 
Ti, the LOFF state has lower free energy than the HS 
state, and the critical behavior of the LOFF momentum 



g ~ (1 - T/Ti 



,1/2 



T 



(39) 



a- 



0.05 



0.04 

ff=0.02 



■••A 



. G 0.8 

T/Ti 



FIG. 4; The LOFF momentum q, scaled by the Fermi mo- 
mentum pf, as a function of T, scaled by the intermediate 
temperature Ti, for several values of number asymmetry a. 



tt=0.02 



0.04 



0.05 



I I 



T/Tf 



FIG. 5: The pairing gap A, scaled by Ep, as a function of 
T, scaled by Tf, for several values of number asymmetry a. 
The solid lines are for LOFF state and the dashed lines for 
BP state. 



In FigUl we show the LOFF momentum q scaled by the 
Fermi momentum pp as a, function of temperature scaled 
by the intermediate temperature Ti for several values of 
number asymmetry a. At any a, the momentum drops 
down from the maximum at T = to zero at T = Ti. 
Note that Ti is a dependent, it increases with increasing 
asymmetry. In Fig[51 we demonstrate the pairing gap A 
for the LOFF state (solid lines) and HS state (dashed 
lines). The LOFF state survives only in the region < 
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T < Ti , and the stable HS state exists in the region Ti < 
T < Tc- The two meet at the intermediate temperature 
Ti with the continuity Aloff(Ti) = /\hs{Ti). Different 
from the pair momentum which increases with increasing 
asymmetry, the gap parameter decreases with increasing 
asymmetry. It is necessary to point out that the LOFF 
state can survive not only in the unstable HS region, but 
also in the region where the HS phase disappears. For 
instance, for a ~ 0.05, the HS phase disappears in the 
region < T < Tq, but LOFF state can survive there. 
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FIG. 6: The free energy scaled by Eq = p|./(87r^m), as 
a function of T, scaled by Tf, for several values of number 
asymmetry a. The solid, dashed and dotted lines correspond, 
respectively to the LOFF, HS and normal states. 



In Fig[B] we compare the free energies for the LOFF, 
HS and normal state for different values of asymmetry a. 
It is clear that LOFF is more stable than the other two in 
the region < T < Ti , and HS is stable than the normal 
state in the region Ti < T < Tc. The phase transi- 
tion from HS to LOFF happens at the intermediate tem- 
perature Ti , and the order parameter characterizing the 
spontaneous rotational symmetry breaking is the LOFF 
momentum q. Since q drops down monotonously to zero 



when T approaches to Ti , the phase transition is of sec- 
ond order. From the above result, the strange superflu- 
idity window in the temperature region To < T < Tc 
obtained infU, disappears, when the LOFF state is 
taken into account. The superfluid can exist in the whole 
temperature region < T < Tc with different phase at 
different temperature. When the temperature is higher 
than the critical value T^ the strong thermal motion sup- 
presses any superfluidity, and the system is in normal 
Fermi gas. From the continuous temperature behavior 
(fTH]) of the HS gap, the phase transition at Tc is of sec- 
ond order. 



ff=0.065 


q/PF 









0.0 02 0.0 04 0.006 0.00 8 0.01 .012 



T/Tf 

FIG. 7: The LOFF pairing gap A, scaled by Ef, and the 
momentum, scaled by pF, as functions of T, scaled by Tf, at 
a = 0.065 > UHS- 





tt=0.065 


LOFF 






Normal \_ 



. 002 0.0 04 0.00 6 0.008 0.01 0.0 12 0.0 14 



T/Tf 

FIG. 8: The free energy scaled by So, as a function of T, 
scaled by Tf, at a = 0.065 > ans- The solid and dotted lines 
are, respectively, for the LOFF and normal phases. 

In the above discussion we have restricted ourself in 
the asymmetry region a < olhs where the HS phase can 
survive. What is about the case with a > a^s, and 
can the LOFF further survive in high asymmetric sys- 
tems where there is no room for the HS? In Fig[7] we 
show the LOFF momentum q and gap A as functions of 
temperature for a = 0.065 > ans- Both drops down 
monotonously with increasing temperature, and the gap 
approaches to zero continuously. To determine the sta- 
ble state, we compare the free energies for the LOFF and 
normal phase in Fig[51 It is easy to see that the LOFF 
is stable at low temperature and the normal phase be- 
comes the only possible state at high temperature, and 
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the phase transition from the normal phase to LOFF is 
of second order. The fact that the LOFF can survive 
at higher density asymmetry where the HS disappears is 
consistent with the conclusion at zero temperature [4^. 
Certainly, when the asymmetry is beyond the maximum 
value aLOFF for the LOFF state, the system is in the 
normal phase at any temperature. 



V. STABILITY OF LOFF PHASE 

In Section IIVI we neglected the possibility of the PS 
phase and focused only on the LOFF phase. Since both 
PS and LOFF phases can be the ground state in the 
temperature region < T < Ti, we analyze in this sec- 
tion the stability of the known LOFF phase against the 
PS. To this end, we study again the imbalance number 
susceptibility x — {d6n/d6^)^ for the LOFF state, 



/d6n\ /9A\ 
fdSn\ / 9g \ 



(40) 



To simplify the expression, we ignore in the following the 
subscript notes and this will not make any confusion. 

Similar to the calculation for the HS phase, from the 
number difference Sn p5|) . we obtain the same expression 
(1^5)1 for the derivatives dSn/dS^ and d6n/dA and the 
derivative 



dSn 
dq 



1 

m 



p • q 

(2^~^ 



[.f\EA) + .nEB)]. (41) 



With the derivatives OA/ 85^ and dq/d6fi evaluated from 
the coupled gap equations for A and q, 



d^n d^n dA d^n dq 



dAdSfi dA^dSfi dAdqdSn 

d^n d^n dA d^n dq 
1 1 . ±. 

dqdS^ dqdA dS^ dq^ dS^ 



= 0, 

= 0, 



(42) 



the other two elements, 

d^n n 1 r dh 

dq^ 



m m? J {2tt) 
+2p-q^ {nEA)-nEB)) 

+q^^{f\EA) + f\EB)) 

p 

1 /■ d^p 



^^:^{f[EA) + f[EB)) 



{l-f{EA)^f{EB)) 



d^fi 



dAdq 771 J (27r) 



pq A 

q Ep 



u'{EA)~r{EB)) 



+q^ {r{EA) + r{EB)) 
p 



+q^{l~f{EA)~f{EB)) 
p 



(45) 



The number susceptibility for the LOFF state is illus- 
trated in Figl9]as a function of temperature in the region 
< T/Ti < 1. For a < 0.53aLOFF, X is always neg- 
ative in the whole region, which means unstable LOFF 
state. This instability of the LOFF state is consistent 
with the conclusion obtained in [l^, 13 ^^^o temper- 
ature. For 0.53aLOFF < a < aHS, X is positive at low 
temperature, then changes sign at the turning tempera- 
ture T2, and keeps to be negative at higher temperature 
between T2 and Ti. Similar to the number susceptibil- 
ity in homogeneous case, it is divergent at the turning 
temperature T2. When a is larger than the maximum 
asymmetry aus for the HS phase but less the maximum 
value aLOFF for the LOFF phase, x is always positive in 
the whole temperature region. 



tt/tlt:=0.28 

ff/tic=0.70 
ff/ffc=0.84 



/I 

' I 

I" 
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II 



we have 



T/Ti 



d6n 
dSfi 



\ OA ) Ikp' 



( 9Sn\' 



I dSn dSn d O 
' OA dq dqdA 



detM 



(43) 

where A4 is the stability matrix for the LOFF state de- 
fined as 



M = 



OA^ dAdq 

d^n d^n 



(44) 



\ dqdA dq^ / 

with the gap susceptibility d^Vl/dA^ given in ([50]) and 



FIG. 9: The number susceptibility x for the LOFF state, 
scaled by its value xo in the symmetric system with a = 0, as 
a function of temperature T, scaled by T\ , for several values 
of number asymmetry a. 



From the analytic and numerical calculations for the 
HS, LOFF, PS and normal state in this section and Sec- 
tions |TT1 [Hi] and |TVl the ground states for the asymmetric 
system at different temperature and asymmetry are sum- 
marized in the phase diagram in T — a plane in Fig llOl 
The dashed line separates the superfluid state at low 
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•.. To Normal 



HS 




a/ttLOFF 

FIG. 10: The phase diagram in T — a plane. The sohd hnes Ti 
and T2 denote, respectively, the first order phase transitions 
from HS to PS and from LOFF to PS, the dashed line means 
the second order phase transition from HS or LOFF to normal 
phase, and F is the tetracritical point. 



temperature and asymmetry from the normal state at 
high temperature and asymmetry. Inside the superfiuid 
state, there are three phases, the homogeneous super- 
fluid above the temperature Ti marked by HS, the stable 
LOFF phase on the right hand side of the temperature 
line T2 marked by LOFF (x > 0), and the unstable LOFF 
phase in between the temperature lines Ti and T2 marked 
by LOFF (x < 0) which is probably the PS phase. Since 
the number susceptibilities for HS and LOFF are both 
divergent at the corresponding turning temperatures Ti 
and T2, the phase transitions from LOFF to possible PS 
and from HS to PS are of first order. On the other hand, 
the phase transitions from HS to normal state and from 
LOFF to normal state are of second order. The four 
phases, HS, PS, LOFF and normal state meet at the 
point F located at (T,a) = {0A9To,aHS = OMaLOFp) 
where Tq is the critical temperature Tc in the symmet- 
ric system with a — 0. The phase diagram we ob- 
tained here for the system with fixed densities is quite 
different from the one for the system with fixed chemical 
potentials ^335 in T — S/i plane. There is no room for 
the PS phase in the T — Sfj, plane, and the window for 
the LOFF and PS states in FiglTU] is much larger than 
the one for the LOFF state in T — 5fi plane. While there 
exists a phase transition line from LOFF to HS in T— Sji 
plane, it disappears in the T — a plane. 

We have checked in the frame of mean field approxima- 
tion that the topological structure of the phase diagram 
in T — a plane does not change for a wide BCS region 
—00 < \/{pFas) < —1- If pair fluctuation is included, the 
second order phase boundary between the normal phase 
and HS or LOFF phase will be shifted, but the topologi- 
cal structure will remain unchanged. Our phase diagram 
is similar to the one obtained in [l3| based on the BdG 
formalism in an atomic trap. 

While our result that the phenomenon of intermedi- 
ate temperature superfiuid shown in Fig[l] is washed out 
by introducing inhomogeneous pairing state is obtained 
by considering the simplest single plane wave LOFF 



state, the qualitative conclusion may remain unchanged 
when we take into account a more complicated LOFF 
state, since a general inhomogeneous state will produce 
a deeper minimum of the free energy of the system. On 
the other hand, we can restudy this problem by inves- 
tigating the stability of the normal phase against an in- 
homogeneous fluctuation [20!. [2ll . Is^ . Let us consider a 
static but inhomogeneous pair fluctuation $(x) for the 
normal state in FigO The effective action of the system 
can be expressed as a series of the fluctuation, 



(46) 



where Sq is the action of the normal Fermi gas. The 
quadratic term S2 reads 



S2 = 



1 



(2^ 



H(k)|$(k)|2, 



(47) 



where $(k) is the Fourier transformation of $(x) and the 
function Ti, is given by 



W(k) = - - 
9 



E 



(27r)3 2e,, -I- |k|2/4m 



(48) 



with Efj^ — e 



p- and = (p ± k/2)^/(2m) - Note 
that Ti. is an even function of k, it depends only on the 
amplitude k = |k|, TL{k). For a symmetric system with 
a = 0, the sufhcient and complete requirement for stable 
normal state is simply expressed as 'H(O) > and the 
critical temperature is determined by 7i(0) = 0, because 
of the condition 7i(fc) > 7i(0). However, this requirement 
becomes incomplete for an asymmetric system with a ^ 
0, since fc = is no longer the minimum of the function 
H(fc). While the equality H(0) = gives two critical 
temperatures To and Tc, and the normal phase below To 
is free from the ordinary BCS instability, characterized by 
7Y(0) > 0, we should demand Ti{k) > for all k to achieve 
a real stable normal state, and the critical temperature 
should be the maximum of those determined by 7i(fc) — 
0. When we calculate the function Ti,{k) in the normal 
phase, for the region close to the critical temperature To, 
7i(fc) is positive at fc = but becomes negative in a range 
of nonzero k. This means that the normal phase there 
is stable against homogeneous superfluidity but unstable 
against inhomogeneous superfluidity. If we consider only 
homogeneous phases, we do find the strange intermediate 
temperature superfluidity. However, this is unrealistic 
when an inhomogeneous condensed phase enters. 



VI. SUMMARY 

We have investigated the temperature behavior of the 
homogeneous and inhomogeneous superfluids in a two- 
component Fermi gas with density imbalance. The main 
conclusions are: 

(l)For homogeneous superfiuid, while the most favored 
pairing temperature is nonzero, and in the case with large 
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asymmetry the superfluidity starts even at finite temper- 
ature, the superfiuid density and the number susceptibil- 
ity are negative at low temperature T < Ti. Therefore, 
the homogeneous superfiuid is stable only at high temper- 
ature T > Ti and unstable at low temperature T < Ti. 

(2) The LOFF phase is energetically more favored than 
the homogeneous superfiuid in the region T < Ti, and 
the phase transition from LOFF to homogeneous super- 
fluid characterized by the pair momentum happens at 
Ti. Due to the existence of the inhomogeneous state at 
low temperature, the so called intermediate temperature 
superfluidityjH, [s^, H, El] disappears and the pairing 
gap A is a monotonous function of the temperature. 

(3) The phase separation is also energetically favored in 
the region where the homogeneous superfiuid is unstable. 
From the calculation of the LOFF number susceptibility, 
the LOFF phase is stable only at large asymmetry, while 
the phase separation should be the ground state at small 
asymmetry. 

(4) When the asymmetry is larger than the maximum 
value for the homogeneous superfiuid and less than the 
limit for the LOFF superfiuid, the homogeneous super- 
fluid disappears at any temperature, but the LOFF su- 
perfiuid survives and is stable at low temperature. When 
the asymmetry is too large, any superfiuid vanishes, and 
the only possible state is the normal state. 



(5)The phase transition from homogeneous superfiuid or 
LOFF superfiuid to normal state is of second order, and 
the transitions from phase separation to homogeneous su- 
perfiuid and to LOFF superfiuid are of first order. The 
four phases meet at the tetracritical point F. The ob- 
tained phase diagram in T — a plane with fixed densities 
is qualitatively different from the one in T — Sfi plane 
with fixed chemical potentials. 

We only analyzed the stability of homogeneous super- 
fluid and LOFF superfiuid against the phase separation 
and neglected the contribution from the surface energy. 
In further works a detailed calculation of the phase sepa- 
ration and the consideration of the surface energvflGllSSj 
are needed. 

Finally we should point out that our result obtained in 
the simple but general two-component model for atomic 
Fermi gas can be extended and applied to other physi- 
cal systems, such as isospin asymmetric nuclear matter 
with neutron-proton pairingfSs', '24', '251 and neutral color 
superconducting quark matter [27, 28]. 
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